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Abstract 

We consider a long Josephson junction excited by a suitable external ac-signal, in order 
to generate control and detect breathers. Studying the nonlinear supratransmission phe¬ 
nomenon in a nonlinear sine-Gordon chain sinusoidally driven, Geniet and Leon m eo 
explored the bifurcation of the energy transmitted into the chain and calculated a threshold 
A(ui) for the external driving signal amplitude, at which the energy flows into the system 
by breathers modes. I numerically study the continuous sine-Gordon model, describing the 
dynamics of the phase difference in a long Josephson junction, in order to deeply investigate 
the “continuous limit” modifications to this threshold. Wherever the energy flows into the 
system due to the nonlinear supratransmission, a peculiar breather localization areas appear 
in a (A, w) parameters space. The emergence of these areas depends on the damping param¬ 
eter value, the bias current, and the waveform of driving external signal. The robustness 
of generated breathers is checked by introducing into the model a thermal noise source to 
mimic the environmental fluctuations. Presented results allows one to consider a cryogenic 
experiment for creation and detection of Josephson breathers. 

The behavior of a nonlinear system excited by an external sinusoidal driving and the bifurca¬ 
tion of the energy transmission into the system were deeply studied by Geniet and Leon mm- 
The Floquet theorem [8j states that, for linear waves in a periodic structure, a forbidden gap 
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of frequency exists. The frequency gap between acoustic and optical branches in the vibration 
modes of a diatomic chain well embodies this picture. Waves with frequency within this forbid¬ 
den range esponentially vanish in the medium (evanescent waves). In nonlinear systems excited 
by plane waves, the appearance of gap solitons [9j and the localization of energy due to non¬ 
linear instability are matter of many systems [25j HOl |23j [13 HU [26]. If the nonlinear medium, 
characterized by a forbidden band gap (FBG), is irradiated to one end by a sinusoidal drive, 
with frequency within the gap, the energy transmission can be supported. This phenomenon is 
called nonlinear supratransmission (NST) and is nowadays reported in several other contexts, 
from Bragg media ED [30], coupled-wave-guide arrays (nonlinear Schrodinger model) [19] 120] . 
optical waves guide arrays [16], Fermi-Pasta-Ulam model ED, Klein-Gordon (KG) electronic 
network [7], chains of coupled oscillators [24] . discrete inductance-capacitance electrical line [18] 
to generic multicomponent nonintegrable nonlinear systems [[]• The manifestation of NST in 
both SG and KG equations, points out that this process is not dependent on the integrability, 
but is a feature of systems possessing a natural forbidden band gap. Geniet and Leon found 
a very simple rule proving an explicit formula for the bifurcation of the energy in a (w, A) di¬ 
agram, called the bifurcation diagram, where uj and A are the frequency and the amplitude of 
the external signal, respectively. The energy flows into the medium as soon as the amplitue A 
of the driving, at frequency w, exceeds the maximum amplitude of the static breather of the 
same frequency. This energy then travels through the system by means of plasma waves and 
nonlinear localized excitations, that is solitons, antisolitons and breathers. Geniet and Leon 
extended the results, obtained for a SG chain, to its continuous version, used to describe a long 
Josephson junction (LJJ) whose extremity is subjected to a sinusoidal excitation ( Neumann 
boundary condition). The soliton solutions are strongly stable, maintaining their shape after 
collisions or reflections, and producing clean evidences in the I-V characteristics of the junction. 
A breather is another travelling SG solution formed by a soliton-antisoliton bounded couple, 
oscillating in an internal frame with a proper internal frequency. The magnetic flux associated 
with a breather is zero, as the voltage difference across the junction. Indeed, considering that 
the time variations of produce a voltage difference across the junction according to the a.c. 
Josephson relation, the rapid oscillations of the breather phase result in a mean voltage close 
to zero (or beyond the sensitivity of nowadays high frequencies oscilloscopes). Other pratical 
difficulties exist: i) a breather has to be efhcently generated and trapped, for a sufRcent long 
time, in a confined area to allow measurements, ii) it decays esponentially in time, even more 
so considering the damping affecting the system, and iii) an applied bias current tends to split 
the breather in a soliton/antisoliton couple. These are the reasons that underlie the lack of 
experimental evidences and confirmations of breathers in LJJ. Observations of breathers, called 
rotobreathers, in Josephson junctions (JJ) are reported only using “Josephson ladders”, i.e. a 
ladder composed by small JJs El El ED ED EH- Our investigation points towards an experimen¬ 
tal setup devoted to the exclusive generation and detection of breathers in LJJ. The influence 
of the environment should be take into account to well describe a real experimental framework. 
Bodo et al E m la enhanced the deterministic analisys developed by Geniet and Leon in a 
discrete SG chain, using a sinusoidal excitation corrupted by an additive white Gaussian noise. 
Considering driving amplitudes below the critical value, they found that noise induced breathers 
can be generated with a probability depending on the noise intensity. 

To check the robustness of the generated breathers, we add to the perturbed SG model a 
Gaussian white source, to mimic the environment effect on the breather’s dynamics. We calcu¬ 
late the percentage of surviving breather, including a white noise source of given amplitude. 

The continuous SG model. — The electrodynamics of a LJJ in the geometrical configuration 
of overlap, is described in terms of a nonlinear partial differential equation for the order parameter 
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<p, that is the sine-Gordon equation mm- Here ip is the phase difference between the wave 
functions describing the superconducting condensate in the two electrodes. Our analysis includes 
a quasiparticle tunneling term and an additional stochastic contribution, ij-]y(x,t), representing 
the Gaussian thermal noise effects. However, the surface resistance of the superconductors is 
neglected. The resulting perturbed SG equation reads 


<p tt (x, t) + Pjip t {x, t ) - tp xx (x, t) = i b (x, t ) - sin(</?(z, t)) + itn(x, *), (!) 

where a simplified notation has been used, with the subscript of ip indicating the partial derivative 
in that variable. This notation will be used throughout all paper. The SG equation is written in 
terms of the dimensionless x and t variables, that are the space and time coordinates normalized 
to the Josephson penetration depth A j and to the inverse of the plasma frequency uj p of the 
junction, respectively. In the limit of small amplitude oscillations, the JJ plasma frequency 
corresponds to the oscillation frequency in the bottom of the a washboard potential well, modified 
by the presence of a bias current [2]. Damping parameter is j3j = (upR^C )^ 1 , where Rn and 
C are the effective normal resistance and capacitance of the junction. The terms i b (x,t ) and 
sin(tp ) of Eq. ([!]) are the bias current and supercurrent, respectively, both normalized to the JJ 
critical current i c . Eq. ([!]) is solved imposing the following boundary conditions 

ip x {0,t) = f(t) ip x (L,t) = 0, (2) 

where f(t) is the external driving pulse used to excite the LJJ and generate the breathers. 

The unperturbed SG equation is 


<Pxx(x, t ) - iptt(x, t ) = sin((^(x, t)). 
and admits solutions in the traveling wave form / = ip(x — ut ) [2], 


ip(x, t) = 4 arctan < exp 




(x — ut) 


(3) 

(4) 


where u is the wave propagation velocity normalized to the speed of light, and is called Swihart 
velocity. Eq. Q represents a single kink , or sohton. The signs + and — indicate a 27r-kink 
(soliton) and a 27r-antikink (antisoliton), respectively. In this framework, ip gives a normal¬ 
ized measure of the magnetic flux through the junction, so that Eq. ([3]) can also represent the 
motion of a single fluxon (or antifluxon). In fact, starting from simple electrodynamic consid¬ 
erations [2], it is possible to obtain a simple relation between the magnetic field H(y ) and the 
spatial derivative of the phase difference 


tp x = 27T 


dH{y) 

<h 0 


integrating over the entire JJ lenght: 


ip(L) — </?(0) = 2ir 


T// 

•ho 


(5) 


where d>o = hc/2e is the fluxon and d = XL+XR+t is the magnetic penetration, A l and Xr are the 
London depths in the left and right superconductors and t is the interlayer thickness. If the phase 
evolution shows a single 27r-kink, a single fluxon will propagate along the junction. Analytic 
expressions of travelling-wave solutions of the unperturbed SG equation can be explicitly written 
also for N-solitons, for kink-antikink collisions and for breathers. Also the damped and biased 
SG equation has slightly perturbed kink as solutions. Indeed, after a series of trasformations in 
zeroth order in bias and damping parameter, the PDE describing the evolution of this perturbed 
system reduces to the unperturbed SG differential equation [IT]. 
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The analytic expression of a stationary breather oscillating with frequency oj is: 


ip(x, t) = 4 arctan 


\/l — oo 2 


sin (oot) 


oj 


cosh ( x\J 1 — oj 2 


( 6 ) 


The SG solutions are invariant under Lorentz transformations, so this stationary breather can 
be boosted into a moving frame, resulting in the following moving breather travelling with an 
envelope velocity v < 1, 


<p(x, t) = 4 arctan 


vr 3 


oj* 


sin [yen [t — vx)] 


OJ 


cosh 


yVl — w 2 (i) (x — 



where 7 = 


(i-»T ,/2 


is the Lorentz factor. 


(7) 


The nonlinear supratransmission: from discrete to continuous SG. — Studying a discrete SG 
chain of N damped oscillators: 

u n — c 2 (u n+ \ — 2u n + u n -i) + sinu n = — r y{n)u n n = l,...,N ( 8 ) 

with boundary conditions 


uo(t) = j4sin(u;£) u n ( 0 ) = 0 u n ( 0 ) = 0 , 


(9) 


where c is the coupling factor and 7 the damping coefRcent, Geniet and Leon [T2f [13] deduced, 
exciting the JJ in the frequencies forbidden gap (w < 1), a threshold .A s (u;): 


A s (oj) = 4 arctan 


—arccosh 
OJ 



1 - oj 2 Y 

UJ )\ ’ 


( 10 ) 


above which NST accours and the system permits the energy trasmission by means of nonlinear 
modes generation (breathers, kinks, antikinks). 

In the continuous limit, that is 1/c —> 0, using the boundary condition: 


d x u\ x=Q = Bsm(oot). 


( 11 ) 


the resulting threshold is: 


B s = 2(1 - oj 2 


( 12 ) 


The B s expression is calculated imposing that the system adapts its phase profile to the breather 
derivative at the boundary d x u\ x=0 = d x Ub [ Varying the position of the breather’s center to 

maximize d x Ub\ x=0 , the simple expression for the maximum amplitude, see Eq. 12 results. 


Numerical results. — The deterministic analysis is carried on setting iTN(x,t) = 0 in the 
Eq. [T| Two different bifurcation diagrams are shown in Fig. [l] Due to the time normalization 
performed for the SG model, set oj < 1 means explore the frequencies lower than the JJ plasma 
one. In the ( A,oj ) parametric space shown in Fig. la the dark spots indicate (A,oj) values 
in correspondence of which only breathers are generated. In particular, the breathers tend to 
localize in branches above the bifurcation threshold, see Eq. [l2j whereas are totally absent 
below the threshold. The grey background of Fig. [la] indicates the absence of breathers modes 
uniquely generated into the junction. In detail, for low (A,lj) values (below the threshold B s ) 
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Figure 1: a) Breathers-branches localized in a ( A , 
medium in a ( A, w) parametric 2D-space. 
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(b) 

') parametric 2D-space. b) Energy injected into the 


there is no NST effect, while above the threshold every kind of nonlinear excitation can propagate 
through the medium. Plasma waves, generated by the external signal and resulting from the 
breathers decay, are also detected with the breathers. Beetwen these breathers-branches, every 
combination of kinks, antikinks, breathers and plasma waves can be detected, but due to the 
stable nature of solitons and the damping of the medium, after a short time only soliton solutions 
survive into the system. The branches for high amplitude values, that is A > 1, corresponde to 
multiple generation of stationary and moving breathers. The results for low frequencies (~ 0.2) 
and high amplitudes are consistent with Geniet and Leon’s results |13j . Changing the bias 
current value, the damping parameter and the duration of the external signal, modifications in 
the breathers arrangement come to light. 

This (A,u) 


The analisys of the energy injected into the system is displayed in Fig. lb 


bifurcation diagram shows the same branched structures, suggesting moreover that the breathers 
tends to localize into isoenergetic regions of this diagram, in particular where the energy is 
slightly lower than the sorrounding areas. Indeed, the energy of a breather [28] : 


-^breather — 16y(l Ld ) 

is less then the energy of a kink/antikink couple 

-S'kink = -S'antikink = 


(13) 


(14) 


The Fig. lb shows also a large white area, sign of the lack of NST (no energy flows into the 
system). 

To check the robustness of the generated breathers, a nonzero noise amplitude is setted. 
The stochastic investigation is developed performing N = 10 3 numerical realizations, setting 
the ( A,u ) values that deterministically corresponde to breathers generation, including into the 
SG model, see Eq. 12 a white noise contribute with amplitude 'Jtn = 10 -3 . The thermal 


fluctuations can supply the system with an amount of energy large enough to create kinks 
instead of a breather, or to break a generated breather into a kink-antikink couple. The results 
in the Fig. [2] shows the percentage of surviving breathers as a function of the amplitude A and 
the frequency lo. For large A values, little amount of thermal energy can be enough to destroy 
a breather. Anyway, a suitable ( A,oj ) values range exists, in which breathers are still generated 
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Figure 2: (Color online) Percentage of surviving breathers as a function of amplitude A and frequency 
u of the external drive, setting 'Ytn = 10 -3 and performing N = 10 3 numerical realizations. 


by external excitation with high probability despite of the thermal influence. 

Conclusions. — We explore the breathers generation and propagation in a long Josephson 
junction externally irradiated by a suitable excitation, changing the amplitude A and the fre¬ 
quency c o of the signal. The analysis is computationally developed in the framework of the 
damped and biased continuous sine-Gordon equation, inserting a Gaussian white noise source to 
include the environmental influence. Taking a cue from the Geniet and Leon results about the 
nonlinear supratransmission in a discrete sine-Gordon chain, we discover a new peculiar localiza¬ 
tion of the breathers in branches in a particular (A, uj) diagram. We distinguish two parameters 
regions: for low values of driving frequencies and amplitudes, no energy flowing into the system 
is detected, otherwise every kind of nonlinear excitations travels through the junction. When 
this occurs, we are able to recognize where only breathers are generated. The deterministic 
analysis is improved including a white noise contribute, to evaluate the percentage of breathers 
enduring the thermal effects. 
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